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Abstract: Let E be the Engel group and D be a rank 2 bracket generating left invariant dis¬ 
tribution with a Lorentzian metric, which is a nondegenerate metric of index 1. In this paper, we 
first study some properties of horizontal curves on E. Second, we prove that time-like normal 
geodesics are locally maximizers in the Engel group, and calculate the explicit expression of 
non-space-like geodesics. 
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1 Introduction 

A sub-Riemannian structure on a manifold M is given by a smoothly varying dis¬ 
tribution D onM and a smoothly varying positively definite metric g on the distribution. 
The triple (M, D, g) is called a sub-Riemannian manifold, which has been applied in 
control theory, quantum physics, C-R geometry and the other areas. Some efforts have 
been made to generalize sub-Riemannian manifold. One of them leads to the following 
question: what kind of geometrical features the mentioned triple will have if we change 
the positively definite metric to an indefinite nondegenerate metric? It is natural to start 

‘Supported by NSFC (No. 11071119, No. 11401531); NSFC-RFBR (No. 11311120055). 

^Corresponding author 


1 



with the Lorentzian metric of index 1. In this case the triple: manifold, distribution and 
Lorentzian metric on the distribution is called a sub-Lorentzian manifold by analogy 
with a Lorentzian manifold. For the details concerning the sub-Lorentzian geometry, 
the reader is referred to [fTSlI . To our knowledge, there are only a few works devoted to 
this subject (see In IfT^ . Chang, Markina, and Vasiliev have 

systematically studied the geodesics in an anti-de Sitter space with a sub-Lorentzian 
metric and a sub-Riemannian metric respectively. In [TTTlI . Grochowski computed reach¬ 
able sets starting from a point in the Heisenberg sub-Lorentzian manifold on K.^. It was 
shown in [12^ that the Heisenberg group H with a Lorentzian metric on possesses the 
uniqueness of Hamiltonian geodesics of time-like or space-like type. 

The Engel group was first named by Cartan [jSlI in 1901. It is a prolongation of a 
three dimensional contact manifold, and is a Goursat manifold. In [[5113IH, A.Ardentov 
and Yu.L.Sachkov computed minimizers on the sub-Riemannian Engel group. In the 
present article, we study the Engel group furnished with a sub-Eorentzian metric. This 
is an interesting example of sub-Eorentzian manifolds, because the Engel group is the 
simplest manifold with nontrivial abnormal extremal trajectories, and the vector distri¬ 
bution of the Engel group is not 2- generating, its growth vector is (2,3,4). We first 
study some properties of horizontal curves in the Engel group. Second, we use the 
Hamiltonian formalism and Pontryagin maximum principle to write the equations for 
geodesics. Eurthermore, we give a complete description of the Hamiltonian geodesics 
in the Engel group. 

Apart from the introduction, this paper contains three sections. Section 2 contains 
some preliminaries as well as definitions of sub-Eorentzian manifolds, the Engel group. 
In Section 3, we study some properties of horizontal curves in the Engel group. In 
Section 4, we prove that the time-like normal geodesics are locally maximal in the 
Engel group , and explicitly calculate the non-space-like Hamiltonian geodesics. 
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2 Preliminaries 


A sub-Lorentzian manifold is a triple (M, D, g), where M is a smooth n-dimensional 
manifold, D is a smooth distribution on M and g is a smoothly varying Lorentzian metric 
on D. For each point p G M, a vector v G Dp is said to be horizontal. An absolutely con¬ 
tinuous curve 7 (t) is said to be horizontal if its derivative y'it) exists almost everywhere 
and lies in Dy(t). 

A vector v G Dp is said to be time-like if g(v, v) < 0; space-like if g(v, v) > 0 or 
V = 0; null(light-like) if g(v, v) = 0 and v 0; and non-space-like if g(v, v) < 0. A curve 
7(0 is said to be time-like if its tangent vector y(t) is time-like a.e.; space-like if 7(0 is 
space-like a.e.; null if 7(0 is null a.e.; non-space-like if 7(0 is non-space-like a.e.. 

By a time orientation of (M, D,g), we mean a continuous time-1 ike vector field on 
M. From now on, we assume that (M, D, g) is time-oriented. If A is a time orientation on 
(M, D, g), then a non-space-like vector v G Dp is said to be future directed if g(v, X(p)) < 
0, and past directed if g(v, X{p)) > 0. Throughout this paper, “f.d.” stands for “future 
directed”, “t.” for “time-like”, and “nspc.” for “non-space-like”. 

Let V, w G D be two non-space-like vectors, we have the following reverse Schwartz 
inequality (see page 144 in [12^ 1: 


|g(v,w)| > ||v|| • ||w|i, 

where |iv|| = yj\g{v, v)|. The equality holds if and only if v and w are linearly dependent. 
We introduce the space Hy(^t) of horizontal nspc. curves: 

Hy(t) = {7 : [0,1] ^ M\ y{f) is absolutely continuous , g{y{f),y{t)) < 0, 

7(0 G Dy(j) for almost all t G [0,1]}. (2.1) 

The sub-Lorentzian length of a horizontal nspc. curve y{t) is defined as follows: 

Ky) = f \\y'{t)\\dt, 

Jo 
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where HyXOII = ^\g{y'it),yit))\- We use the length to define the sub-Lorentzian dis- 
tanee duiquqi) with respeet to a set t/ c M between two points q\,q 2 e U: 

f sup{liy),y e Hu(qi,q 2 )} if Hyiquqi) * ® 
du{qi,qi) = \ 

0 if Hu{qi,q 2 ) = fb, 

where Hu{qi,q 2 ) is the set of all nspc.f.d curves contained in U and joining q^ and q 2 . 

A nspc. curve is said to be a maximizer if it realizes the distance between its end¬ 
points. We also use the name t/-geodesic for a curve in U whose each suitably short 
sub-arc is a [/-maximizer. 

A distribution D c TM is called bracket generating if any local frame {X,}i<,<r for 
D, together with all of its iterated Lie brackets [X,, [Xj,Xk]], • • • span the tan¬ 

gent bundle TM. Bracket generating distributions are sometimes also called completely 
nonholonomic distributions, or distributions satisfying Hdrmander’s condition. 

Theorem 2.1. (Chow) Fix a point q £ M. If the distribution D c TM is bracket 
generating then the set of points that can be connected to q by a horizontal curve is the 
component of M containing q. 

By Chow’s Theorem, we know that if D is bracket generating and M is connected, 
then any two points of M can be joined by a horizontal curve. 

Now, we describe the Engel group E. We consider the Engel group E with coordi¬ 
nates q = {x\,X 2 ,y, z) 6 The group law is denoted by O and defined as follows: 

ixi,X2,y,z)Q(xfx2,y,z!) 

, , , XiX2-x\X2 , X 2 X 2 , XiX' 

Xi + Xi,X 2 + X 2 ,y + y + - - -,z-i-z + ^-(x 2 + X 2 ) + xiy + ^-(xi + x^) 

A vector field X is said to be left-invariant if it satisfies dLqX{e) = X(q), where Lq 
denotes the left translation p Lq(p) = qQp and e is the identity of E. This definition 
implies that any left-invariant vector field on E is a linear combination of the following 
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vector fields: 


^ d X 2 d ^ d ^ x\ d ^ x\ + x\ d 

^ dx\ 2 dy" ^ dx 2 2 dy 2 dz 

d d d 

oy oz oz 


( 2 . 2 ) 


The distribution D = span{Xi,X 2 } of E satisfies the bracket generating condition, 
since X 3 = \Xx,X 2 \,Xn = [Xi,X 3 ]. The Engel group is a nilpotent Lie group, since 
[Xi,X 4 ] = [^ 2 ,^ 3 ] = [^ 2 ,^ 4 ] = 0. We define a smooth Lorentzian metric g on E such 
that g(Xi,Xj) = (-\y''6ij, i,j = 1, - • ■ ,4, where d,y is the Kronecker symbol. It is not 
difficult to compute the coefficients of g under the local coordinates (x\,X 2 ,y,z) 6 R'*. 
The coefficients can be expressed as 
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(2.3) 


When we restrict g to D, we can get a smooth sub-Lorentzian metric g 
satisfies 


go, which 


g(Xi,Xi) = -1, g(X2,X2) = 1, g(Xi,X2) = 0. (2.4) 

On the other hand, any sub-Lorentzian metric on D can be extended to a (usually not 
unique) Lorentzian metric on E. In this paper, we assume that Xi is the time orientation. 

3 Horizontal curves 

Chow’s theorem states that any two points can be connected by a horizontal curve, 
but we have no information about the character of horizontal curves. In this section, we 
will investigate some properties of horizontal curves. 
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An absolutely continuous curve 7 ( 5 ) : [0,1] ^ £ is said to be horizontal if the 
tangent vector 7 ( 5 ) can be expressed linearly by the horizontal directions X\, X 2 , hence 
we have the following lemma. 


Lemma 3.1. A curve y(s) = (xi(s), X 2 (s),y(s),z(s)) is horizontal with respect to the 
distribution D, if and only if 


X2X\ 

~Y' 


X\X2 

~ir 


+7 = 0, 


x1 + X 

~ 


2 

—X 2 + z = 0. 


(3.1) 


Proof The distribution D is the annihilator of the one-forms: 

X2 Xi Xj -I- x\ 

oji = —dx\ - —dx2 + dy, <x>2 =-^—^'*‘•2 + 

so 7 ( 5 ) is horizontal if and only if (13.11) holds. 


□ 


By the same method, we can easily calculate the left invariant coordinates uds) 
and U 2 {s) of the horizontal curve y(s): 

U{ = x\, U 2 = X2. (3.2) 

The square of the velocity vector for the horizontal curve is: 

g{y, y) = -u] + ul = -Vi + i 2 - (3.3) 

So whether a horizontal curve is time-like(or nspc.) is determined by the sign of -x\+x\. 

Next we present a left invariant property of horizontal curves in the Lie group 
with sub-Lorentzian metric. That is to say, the causal character (time-like, space-like, 
light-like, or non-space-like) of horizontal curves will not change under left translations. 
Hence it is also true for the Engel group. 

Let us consider a left-invariant sub-Lorentzian structure on a Lie group G: ^ = 
span{Xi,X 2 ,-■ ■ ,Xk) TG, g{Xi,Xj) = (-l)^‘'d,j, with a time orientation Xi. The 
vector fields X, are left-invariant, i.e. 

Lx*Xi(q) = Xfx ■ q), x,q £ G, i = I, - ■ ■ ,k. 
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Proposition 3.2. Left translations preserve the causal character of horizontal curves of 
a left-invariant sub-Lorentzian structure on a Lie group G, and the property of future¬ 
directness is also preserved. 

Proof Let c{f) be a causal horizontal curve, and 

k 

c(t) = ^ Ui(t)Xi(cit)). 

i=l 

Then, the left translation yit) = xO c(t) has the same causal character, since 

k k 

y(t) = L„c(t) = LUJ] Ui(t)Xfcm = Yj Ui(t)LMicm 
1=1 1=1 

k k 

= ^ UiiOXfx O c(0) = ^ Ui(t)Xi(y(t)). 

1=1 (=1 

Therefore, 

k 

g{c{t),c{t)) = ^(-l)‘^'‘Mf = g{y{t),y{t)), 

1=1 

g(c(0,^i) = -Ml = ^(7(0, ^l)- 


□ 


By Chow’s Theorem, we know that any two points on the Engel group can be 
connected by a horizontal curve. But we do not know its causal character(time-likeness, 
space-likeness, light-likeness). This is not an easy problem. We are able to present some 
particular examples to show its complexity. 

Example 1: Let X 2 = 0. Then X 2 = x\ is constant. The horizontal condition (13.11) 
becomes 


X2 

y Xi +7 = 0, 

(3.4) 

z = 0. 

(3.5) 
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And the square of the veloeity vector 

-u\ + u\ = -x\ < 0. (3.6) 

It follows that, the curves satisfying (13.41) and (13.51) are all non-space-like curves. Fur¬ 
thermore, we obtain, 

y(5) = -^X 2 Xi(s) -H ^X 2 X^^ + y°, z(5) = z°. (3.7) 

Therefore, all nonconstant horizontal curves c(s) = (vi(5), x^, -^ + y°, z°) are 

time-like. These curves are straight lines. If ii = 0, c(5') degenerate into some points, 
so there are no null curves in this family. 

Example 2: Let X 2 0. We choose X 2 as a parameter, then the horizontal condition 
(13.11) becomes 


•^2 . .^1 . „ 

(3.8) 

x] + xl 

- 2 = 

(3.9) 

And the square of the velocity vector 


—Wj + U2 — —+ 1. 

(3.10) 

We consider there different cases. 


(a) If x\ = 0, then xi = x^ is constant, (|3.8|) and (l3.9l) become 


-^+^ = 0. 

(3.11) 

- 2 = 

(3.12) 

In this case. Icf.vfp = 1. so the curves satisfvine (13.111) and (13.121) 

Furthermore, we obtain, 

x^ 1 (''• 1 ) 

+ /, Z(S) = ,4 + 0 ■>^2 + Zo- 

1 0 2 

are all space-like. 

(3.13) 
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Therefore, all nonconstant horizontal curves c(s) = (x^, X 2 , + ^^-^X 2 + Zo) 

are space-like. There are no null or time-like horizontal curves in this family. 

(b) If y = 0, (13.81) and (13.91) become 


X2X1 - xi = 0 , 


-I- xl 


-I- z = 0. 


(3.14) 

(3.15) 


From (I3.14L we get 


1 

JC2 




integrating with respect to X 2 , we calculate Xi 
substituting xi in (13.151) . we obtain 


lx 2 , where l = -i, i.e. Xi = -ix 2 . 




xl + z°. 


Therefore, all nonconstant horizontal curves 


c(5) = \lx2 


>,.^ 2 ,/, ^(1 +L^)xl + z^'^ 


(3.16) 


(3.17) 


are time-like when |t| > 1. If |t| < 1(= 1), they are space-like(null). 
(c) If z = 0, the horizontal condition becomes: 

^2 . „ 

Y,-j+y = o, 


X j -I- xl 


0 . 


(3.18) 

(3.19) 


So .Vi = .^2 = 0, y = yo- The curves degenerate into some points. There are no causal 
(time-like, space-like, null) horizontal curves in this family. 

Thus, any two points /’i(A:'j’, V 2 ,y°,z°), Qi(xi,xl,y^,z^) can be connected by a time- 

XiX^ ^ 0^0 

like horizontal curve if + y^. Especially, any two points (Xp0,y^,z^), 

(xi , 0, y°, z°) can be connected by a time-like horizontal straight line. 

Any two points xl,y^,z^), X 2 ,y'^,z^) can be connected by a space-like 


horizontal curve if y' = yX 2 + y“, z' = \x\ + + Zo- 


(-?)^ 
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Any two points Piix^, y°, z°), Q^ixi , can be connected by a time-like(space- 


,0 


like, null) horizontal curve if x\ = 1 x 2 , ^ (l + xl+z°, and \ l\= ^ > 1(< 1, = 1). 

4 Sub-Lorentzian geodesics 


In the Lorentzian geometry there are no curves of minimal length because two 
arbitrary points can be connected by a piecewise light-like curve whose length is always 
0. For example, let be the two dimensional Minkowski space, p = (x, y) is any one 
point in this space. We want to find a light-like curve going from the origin to p. First, 
we choose a curve 71 (t) : (x(t),y(t)) = (t,t) which connects the origin and the point 
(£|y, choose the second curve 72(0 : {x{t),y{t)) = {t,-t + x + y) which 

joints (^, and p. It is easy to check that the curve y{t) consisting of 71 and 72 is 
a light-like curve. It goes from the origin to the point p, and the length is 0. However, 
there do exist time-like curves with maximal length which are time-like geodesics 
Upon this reason, we will study the optimality of time-like geodesics, and compute the 
longest curve among all horizontal time-like ones on the sub-Lorentzian Engel group. 
The computation will be given by extremizing the action integral S = ^ f (-u^ -l- u\)dt 
under constraint (13.11) . By Proposition 13.21 horizontal time-like curves are left invariant, 
so we can assume that the initial point is origin, i.e., Xi( 0 ) = X 2 ( 0 ) = 7 ( 0 ) = z( 0 ) = 0 , 
and time-like initial velocity is -u\{0) + u\{0) = - 1 . 

Let ^ = (^ 1 , ^ 2 , ^ 3 , ^ 4 ) be the vector of costate variables, so the Hamiltonian function 
of Pontryagin’s maximum principle is 

X\U 2 - X2U1 


2 2 

= ^0 ——- + ^\U\ +^2^2 + ^3- 


x^ -I- xl 

+ ^4- " 


U 2 . 


2 2 ""2 

where ^0 is a constant equals to 0 or -1. Also, we get the Hamiltonian system: 

2 2 

X^lli — XoU] -^1 -^9 

Xl = = Uu X2 = //^2 = ^ 2 , = ^^3 = -^^ Z = = —-— 


(4.1) 


U2, 


= -H. 


6 “2 


■^1 


-^4X1^2, ^2 = - H _ 


6m 1 


^2 


- 6-^2M2, 6 = 6 = 


(4.2) 
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and the maximum condition: 


H{^q, q{t), u{t)) = max H{^q, q(t), u), < 0, (4.3) 

U€R2 

where u(t) is the optimal control, and ^(0) ^ 0- 


4.1 Abnormal extremal trajectories 


We shall investigate the abnormal case = 0- From the maximum condition (14.31) 
we obtain 


= 0 , 


H. 


U2 


6 + 


6^1 ^4(^1 + x\) 


+ 


0 . 


2 2 

Differentiating equations (14.41) and (14.51) . we obtain 

0 = ^ 1 -%^=^!-^ = + ^4^1), 


0 = 6 + 


6-^1 


+ ^ 4 (XiXi + X2X2) - Mi (6 + ^4X1). 


(4.4) 

(4.5) 


(4.6) 

(4.7) 


For the time-1 ike curve, we assume that -Mj -l- ^2 = -1, so 6 + 6-^i = 0. If 6 = 0, 
then 6 = 0, and therefore ^ = 0. It is a contradiction with the nontriviality of the costate 
variables, hence 6 ^ 0. In this case, Vi = ^ is a constant, and Ui =0, U2 = +i, so 
there is no time-like abnormal extremal in the Engel group E. 

For the space-like curve, we assume that -Mj -i- = 1, by using the same method, 

we get that ui =0, M 2 = ±1, so the space-like abnormal extremal trajectories are given 
by the following expression: 


y(s) = 


0, +5', 0,+—|. 


(4.8) 


For the null curve, suppose that -Mj -l- =0, we can easily get that ui = 0, U 2 = 0, 

so the null abnormal extremal trajectories are trivial curves. 
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4.2 Normal geodesics 

4.2.1 Normal Hamiltonian system 


Now we look at the normal ease = -1 • It follows from the maximum eondition 
(14.31) that Hui = Hu 2 = 0. Henee 


Ui = - 




^3x1 ^4(^1 + 

^2 — s2 —;;—I-- 


(4.9) 


Let = (^, Xj),i = 1,2, be the Hamiltonian eorresponding to the basis veetor fields 
Xi,X 2 in the eotangent spaee T*E. They are linear on the fibers of the eotangent spaee 
T*E, and 

X? Xi X? + xi 

6 = 6 + y6 + (4.10) 

So Ml = -^1 and M 2 = 6- 

The Hamiltonian system in the normal ease beeomes: 

' i, = f = -(6 - f 6 ) = - 6 , 

6 = = (6 + §-6 + 

“ if “ = 2(-^*6 + -^26)5 


- “If - “6(| + -^16), 
6 = -If = -^66 - ^266, 

ft = -f = 0, 

.^4 = -f = 0. 


Definition 4.1. A normal geodesic in the sub-Lorentzian manifold (E, D, g) is a curve 
y : [a, b] ^ E that admits a lift Y : [a, h] ^ T*M, which is a solution of the Hamiltonian 
equations ( l4.iiD . In this case, we say that Y is a normal lift ofy. 


Associate with the expression of H, a sub-Lorentzian geodesic is time-like if // < 
0; space-like if // > 0; light-like if // = 0. 
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Remark 4.1. In fact, abnormal extremal trajectories d4.^l) are also normal geodesics, 
since we can choose the costate variables as ^ = (0, +1,0,0), it is easy to check that 
r(0 = (Tj^) satisfies Hamiltonian equation M.llk This example also confirms that 
normal geodesics and abnormal trajectories are sometimes not mutually exclusive. 

Lemma 4.2. The causal character of normal sub-Lorentzian geodesics does not depend 
on time. 

Proof. The Hamiltonian H is an integral of the Hamiltonian system, i.e., H{s) = 0, this 
implies that the eausality eharacter does not change for all t 6 [0, oo). □ 

Remark 4.3. Ify{t) is a nspc. normal geodesic on the Engel group, then the orientation 
will not change along the curve. In fact, ifyit) is time-like, and it is future directed at 
t = 0, then we have -u\{t) + u\{t) = -1, ufQ) > 0. We only need to show that uft) will 
not equal to 0 along the curve yit). Actually, if there is a t\ > 0, such that mi(0) = 0, 
then we have u\{ti) = -1, it is impossible. So u\{t) will not change the sign (since 
uft) = -^\{t) is continuous), and y(t) is future directed along the curve. It is also true 
for the other cases. 

4.2.2 Maximality of short arcs of geodesics 

Definition 4.2. Let (pbe a smooth function on M, U is an open subset in M, the horizon¬ 
tal gradients! of(p is a smooth horizontal vector field on U such that for each p £ U 
and v£H, d„(p(p) = gC^H^Pip), v). 

Locally, we can write 

r 

^ hT = -idxiP^)Xi + fdxjfi)Xi. 

i=2 

Now we give a proof that the time-like normal geodesics are locally maximizing 
curves on the Engel group. 
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Proposition 4.4. Ify is a t.f.d. (t.p.d.) normal geodesic on the Engel group, then every 
sufficiently short subarc ofy is a maximizer. 

Proof. Assume that y : (a, b) ^ E is parameterized by are-length, y(t) = u^ft)Xi(y(t)) + 
U 2 (t)X 2 (y(t)), Xi is the time orientation, and f (t) = (y(t), A(t)) is the normal lift of y. 
So we have H(y(t),A(t)) = t 6 (a,b). For any c 6 {a,b),e > 0, let Jc = (c - 
6,c + e) c (a,b)he a neighborhood of c. We will prove that y\j^ is maximal for any 
c e (a, b) and small 6 > 0. Sinee the sub-Lorentzian metrie is left invariant, so we 
can assume that y(c) = 0, A(c) = Aq. Consider an n - 1 dimensional hypersurface S 
passing through the origin 0, and satisfying /lo(7’o(5')) = 0. Let d be a smooth one-form 
on an open neighborhood Q of 0, such that d(0) = do, and Vp e S n Q, A(p)(TpS) = 0, 
H{p,A{p)) = -i. Let Yp = (yp,Ap) be the solution of r(0 = f^(r(t)), r(c) = (p,d(p)). 
Then clearly Fq = f. Since 7(0) ^ TqS, by the Implicit Function Theorem, there exits a 
diffeomorphism: 

V : {c - e,c + e)xW ^ U <z E, 

{t, p) yp{t), 

where VF is a neighborhood of 0 in S, t/ c Q is a neighborhood of 0 in Define a 
smooth function V : U R as: 


V(x) = t, if X = yp(t), 

we will show that ||V//y|| = 1. For this purpose, let be the vector field on U given by 

= JpiA) = ufp, fXfypit)) + U 2 {p, t)X 2 {yp{t)), if x = yp(t), 

where ufp, t), U 2 (p, t) are smooth functions on VFx(c-e, c-l-e), and Mi(0 , t) = u^ft), U 2 ( 0 , t) 
u^it). Since H(p,A(p)) = by the construction of Fp(0, we have//( 7 p(t), dp(0) = - 5 , 
and -Mj + ul = -1. It is easy to check that Ti = uiXi + F 2 = U 2 X 1 + U 1 X 2 
is also an orthonormal basis of D, so = 1, dy^^V = 0. Therefore, VhV = -Ti, 
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||V//y|| = -^01 = -^1 - Mj + u\\ = 1. Choose t\,t 2 in the domain of y. If y{t) is 

a t.f.d. geodesic, then \u ?l > 1 ^^°!, and Mj > 0. Since Mi(0, t) = Wj, and Ui(p, t) is a smooth 
function, so there exists a neighborhood Wi x (c - 6 i, c + 6 i) c Wx{c-e,c + e) such that 
u\{p, t) > 0. Thus VhV = -Y\ is past directed. On the other hand, since -u\ + u\ = 
we have \u\\ > \u 2 \. Let 77 : [0, a] ^ [/ be a t.f.d. curve with ? 7 ( 0 ) = y(ti), 77 ( 0 ;) = y{t 2 ). 


and ?7 = ViXi + V 2 X 2 , then |vi| > IV 2 I, Vi > 0 , so g{f],X hV) = UiV\ - U 2 V 2 > 0 , and 

dV(r](s)) 


^rb,, 7 .]) = t 2 -h = V(y(t 2 )) - V(y(h)) 


-f 


lb 

0 


/ 


ds 


-ds 


g(.f],VHV)ds > \\f]{s)\\ds = L{p\[o^a])- 


By the reverse Schwartz inequality, L(y) = L{7]) holds if and only if 77 can be repa¬ 
rameterized as a trajectory of -VhV. If y(t) is a t.p.d. geodesic, then |i 7 j| > b^l, and 
< 0. By the same method, we choose a neighborhood such that 1^2 x (c - 62 , c - 1 - 62 ) c 
W X (c - e, c -I- e) such that ui(p,t) < 0. Thus VhV = -Ti is future directed. Let 
p : [0, or] ^ t/ be a t.p.d. curve with p(0) = y(ti), p(a) = y{t 2 ), and p = piXi - 1 - P 2 X 2 , 
then Ipil > |p 2 l,yni < 0, so g(p, VhV) = WiPi - U 2 P 2 > 0, and 

dVijiis)) 


= t 2 -h = V{y{t 2 )) - V{yih)) 


J r*0! r*Q 

g{p,VHV)ds> 

0 Jo 


-f 


ds 


-ds 


Iip(5)||j5 = L(p|[o,ff]). 


By the reverse Schwartz inequality, L(y) = L(p) holds if and only if p can be reparam¬ 
eterized as a trajectory of -VhV. In conclusion, the t.f.d(t.p.d.) normal geodesics are 
locally maximizers. This ends the proof. □ 


Next, we compute the expressions of light-like geodesics and time-like geodesics 
on the Engel group. 

Differentiating 


^1 - ^1 “ y-^2 - -^ 2(^3 + -^ 1 ^ 4 ), (4.12) 

^2 = ^2 + + {X\Xi + X2i2)^4 = + •^1^4)- (4.13) 
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Let 


A^) = -(^3 + ^i^4), li = Uu (4.14) 

then we have 

L =K2. 6=ye^i, $ = (4.15) 


4.2.3 Light-like geodesics 

Firstly, we study the ease of light-like sub-Lorentzian geodesies. 

By the definition, we have H = \{-h\ -l- h^) = 0, thus /z 2 = ±fii. If ^2 = ^ 1 , then 
light-like trajectories satisfy the ODE: 

y = -h,{X,-X2), 

i.e. they are reparameterizations of the one-parametric subgroup of the field X\ - X 2 . 
We assume j = Xi - X 2 , so 

1 1 2 2 

ii = I , X 2 = -1, y = --(21:1 + X 2 ), z = --(21:1 + 21:2)’ 

thus 

^ 1 3 

Vi = t, X 2 = -t, y = 0, z = --r. 

If h 2 = -hi, similarly, we obtain 

^ 1 3 

xi = t, X 2 = t, y = 0, z = -t . 

In conclusion, we get the following theorem: 

Theorem 4.5. Light-like horizontal geodesics starting from the origin are reparameter¬ 
izations of the curves: 

1 3 

.vi =t, X 2 = +t, y = Q, z = ±-r, 

i.e., they are reparameterizations of the one-parameter subgroups corresponding to the 
vector fields X\ + X 2 . 
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4.2.4 Time-like geodesics 


Secondly, we study time-like sub-Lorentzian geodesics on the Engel group. 

We consider the case of ^4 = 0 at first. This case is also of interest since it re¬ 
produces the earlier known results for the Heisenberg group [[2^ . In this case = 
-(^3 -I- = -^3 is a constant. Equations (14.151) become 

6 = - 6 ^ 1 , (4.16) 


where ^3 is a constant. There are two separate cases: 

Case 1: If ^3 = 0, we have and ^2 are constants, i.e., ^ 1 ( 5 ) = ^i(O) = ^i(O) and 
= ^ 2 ( 0 ) = ^ 2 ( 0 ). According to (??), and ^2 are constants. On the other hand, by 
integrating ii = -^1 and X 2 = ^ 2 , we get 


Vi( 5 ) = -^is and X 2 (s) = ^ 2 ^- 


(4.17) 


Since y = kxi ^2 + .*^ 2 ^ 1 ) = 0, then y(s) = 0. Also 


z = 


-|- y^ ^ 4- ^ 

■^1 ^ "^2 y Si ^ S2 . 2 

——6 = —, 


so 

, , 3 x]{s)X2{s) + xl(s) 

z(.) = —ft. =-^-. 

Theorem 4.6. In the case of ^2 =^4 = 0. there is a unique time-like horizontal geodesic 
joining the origin to a point {x\,X 2 ,y,z), if and only ify = 0, z is the following function 
ofxi,X2: 


X^,X2 -I- xl 

z = —^- - 

6 

And the expression of the geodesic is 


(4.18) 


Xi(5) = -^i5, X 2 {s) = ^ 2 S, y{s) = 0, z(s) = 



(4.19) 


17 










where , ^2 are constants. The arc-length is given by 


I = - x^. (4.20) 

Its projection to the {xi,X 2 ) plane is a straight line. 

Case 2: If ^3 4^ 0, from (14.161) . we have 

^ 1 ( 5 ) = eosh(^35) - sinh(^35), (4.21) 

= -t\ sinh(^35) + eosh(^35), (4.22) 


where = ^i(O), = ^ 2 ( 0 ). So 

Vl 


iO ^ 

- I ljx(t)dt = -^ sinh(^35) + (eosh(^35) - 1), 

Jo 


ft 


ft 


V2 


ft ft 

I ft( 0 ^i?? = - 7 ^ (cosh(^35) - 1) + ^ sinh(^35). 
Jo 


(4.23) 

(4.24) 


/o ft ft 

Substituting them into the expression of j,z in (I4.11L and integrating, we get 

Theorem 4.7. In the case of^j ^ O^ft = 0. fhe time-like horizontal geodesics starting 
from the origin are given by: 


xi(5) = -Ai sinh(^35) + A 2 (eosh(^35) - 1), (4.25) 

X2(s) = -Ai (eosh(^35) - 1) + A 2 sinh(^35), (4.26) 

y(s) = ^(Al - A?) (^3^ - sinh(6^)), (4.27) 

zis) = A 2 (Ai + aI) eosh^(^35) sinh(^35) - ^aI sinh^(^35) - ^^ 1(^1 + cosh^(^35) 

+ ^Ai(AI + 3^2) cosh'(^35) - ^A2(3 Ai + A2) sinh(^35) cosh(^35) - ^^2(3^^ -v aI)s 

-^AdAl + 3Al). (4.28) 

o 

where = ^i(0),^2 = ft(0) is the initial value, ^3, Ai = ^, A2 = ^ are constants. 

Projections of geodesics to the plane (xi , X2) are hyperbolas, for ^( 0 ) = ( ^^ 2 , 1 , 1 , 0 ), 
^( 0 ) = (-^, 1 , 0 ) and ^( 0 ) = (-^, 5, - 1 , 0 ), they are shown in Figure[IJ 
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Figure 1: Projections of geodesics to the plane (xi,X 2 ) when ^3 0, ^4 = 0 


From this theorem, we obtain a description of the reachable set by geodesics ^3 3 ^ 
0 , ^4 = 0 starting from the origin. 

Corollary 4.8. In the case of ^3 3 ^ 0,^4 = 0, let (xi,X 2 ,y,z) be a point on a time-like 
geodesic, then we have 


-1 < 


4y 


-Xj + x\ 


< 1 . 


Proof. By (14.251) and (I4.26L we get 


f-i 


-Xj + ^2 = 4(^2 - Aj) sinh j, 

substituting (14.291) into (I4.27L we obtain the following equation: 

{-x] + xl ){,^3 - sinh(^ 3 )) 


3 ^ = 


8 sinh^ (I) 


(4.29) 


(4.30) 


if we set T = y, then 


or 


{-x\ + x]) 


Ay 


sinm(T) 


- coth(T) , 


- COth(T). 


(4.31) 


(4.32) 


(-Xj + x^ sinh^(T) 

It is easy to check that the right hand side of (14.321) is a decreasing function in (- 00 , + 00 ), 
and its range is (-1,1). That is to say, the points on the time-like geodesics should satisfy 

4y 


-1 < 


-I- x\ 


< 1 . 


This ends the proof. 


□ 


19 































Next, we consider the case ^4 0 . Recall that 


=P^ 2 , 6 =^^i, where ^( 5 ) = -(^3 + ^ 1 ^ 4 ), $ = ^ 4 ^i- 
Combining the expressions for yS and ^2 to get 


U2 =/3U^i =m 


Integrating both sides, we have 


«2 

^46 = — + Cl, where Ci = ^ 4 ^( 0 )-z— = ^ 4^2 “ 


This yields 


and 


Xi(5) = 


>e(^) + 6 

^4 ^ 




^4\ 2 


Since X 2 = ^ 2 , we deduce 

X2(5) = f ^2it)dt = ^ f 

Jo S4 Jo 


■s /a2 




+ Cl I dt. 


To compute yCi') in term of /3{s), we note that 

y = ^iXl ^2 + X 2 ^l) = ^ixiX2 - X2X1), 


then integration by parts yields 

1 r 

y(^) = 2 J {X1X2 - X2Xi)dt 


-f 


X\^2dt - -X1X2 


m)+a)(^^ + Cijdt-^xiX2. 


Finally, since i; = ^^^ 2 , 


+ 1^9 1 

" J — " 2 Jo ^ 6 


xX 






(4.33) 


(4.34) 


(4.35) 


(4.36) 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


(4.41) 
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Once we find jS, we can find the geodesic (a:i(5), X 2 {s), y{s), explicitly. 
Sincey8(5) = ^4^i,y8(0) = = ^ 4 ^°’ we have 

/3\s) 


m = ^ 4^1 = =m{U2) =m 

Multiplying both sides by 20(s) and integrating, we have 


+ Ci . 


$\s) = ^ + Crfi^s) + C2 = + Cl I + C2 - Cl 


where C 2 is a eonstant, and 


C 2 = 0lO) - ^ - CyfilO) = ieiff, + I - 


Then 


0 


C2-Ct = (OY 4 + J- ^2^3^4 - I ^2"^4 - ^ I = ^4((^1 )" " (^ 2 )') = 0 

sinee (^)2 - (^0 = 1. 

Assume 0(s) > 0, we have 


^3 


^2//^0x2 //r0\2x 


54’ 


dfijs) ^ i(0\s) 

ds 


+cii +e 4 - 


Henee 


ds 


d/3 


V{¥^ 




(4.42) 


(4.43) 


(4.44) 


(4.45) 


(4.46) 


(4.47) 


Let = Cl + ^ 4 /, p^ = Cl - ^ 4 / and u = integrating (14.471) from 0 to s, we obtain 


= r 

Jm 


PU) 

V2 


y/ldi 


Set 


^2 ^ (p - p) 


^ V(m2 +p2)(M2 +p2) 

/cl+Pl-C, 


(4.48) 


4pp 


■ V^;+4 ’ 


g = 


2V^ 2(C2 + ^2)i 
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Since 


r 


dt 




= g • cn ' {cos ip, k) = g,F{p, k). 


(4.49) 


where cn ' (jy, k) is a Jaeobi’s Inverse Elliptie Funetions, and 


^ = eos-i (4-^)’ F{p,k)= r 
y^+PP) Jo 


dt 


Vl - k^sird t 


Henee 


X 


«£) 

VI 


VIjii 


m 

VI 


+ p^)(ld + p2) 
Aeeording to (14.491) . we have 


= r 


V2<ij( 


AO) 

VI 


^OO 

JM 

VI 


V2Jii 


V(m^ +p2)(M2 +p2) Jm y/(M2 +p^){u^ +p2) 


where 


Sinee 


Henee 


f 

J&. 


y/2du 


m 

VI 




= V2gF(^i,fc) = constant. 


pi = cos 




^^2dl 




= V2g- 


cn 


-1 (^is)j-2£p 
yS2(5) + 2pp 


-1 ifis)-2pp\ 5 


let F = F{pi,k), we obtain 


(4.50) 


(4.51) 


(4.52) 


ySX^) = 


2 pp(1 + cn (F - 2pp(1 + cn (25, k)) 

- cn {2~s,k)) ’ 




(4.53) 


where 2s = F — 
VIg 
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Since 


\ - cn ( 2 s) 
\ + cn ( 2 s) 


tn^ (s)dn^ (s). 


(4.54) 


hence 


m = 


sj2^ 


= sj 2 ppcs (S, k)nd (s, k). 


(4.55) 


tn (s, k)dn (s, k) 

For the case of 0(s) < 0, we can calculate by the same method, and get the same 
result. But the expression of the parameter s in (14.531) and (l4.55l) should be changed to 


^|F + 


VIg 


Thus the sign of $(s) will not affect the expression of the geodesics. 

Therefore, integrating equations (14.361) . (I4.38L (14.401) and (14.411) . we get a complete 
description of the Hamiltonian time-like geodesics in the Engel group. 


Theorem 4.9. In the case of ^4 ^ 0, time-like geodesics starting from the origin are 
given by: 


.^i('^) = --^(^(5) + ^ 3 ), 

.^2(^) = :^(B2(s) + 2Ci), 

2^4 


1 


y(s) = -:^(B3 (s) -V 2 CiBi(s) -h ^3^2(5) + 2 Ci^ 35 ) - ^Xi(s)x2(s), 


(4.56) 

(4.57) 

(4.58) 


1 


1 


Z(s) = —AB 4 (s) + 2 C,B 2 (s) + 2 ^ 353 ( 5 ) + 4Ci6Bi(5) + fiB 2 (s) + 2Ci^3^) + 7 ^ 2 (^). 


(4.59) 


where C\ = ^ 4 ^° “ "f- Bfs) = f3‘(t)dt, i = 1,... ,4, and the expressions of Bfs) are 

presented in Appendix. 


Projections of geodesics to the plane (xi,X 2 ) with ^(0) = (1,0,1,1), ^(0) = 
(-^, 1/2,2,1) and ^(0) = (-^, 1/2,1,1) are shown in Figure [2l 
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Figure 2: Projections of geodesics to the plane {xi,X2) when ^3 0, ^4 0 

Appendix 

Denoting 


1/^1 ( 5 ) = In (k'^ + cs^ (s, k)). 


if'iis) = 


dn {2s, k)sn {2s, k) 


^^{s) 


E{2s, k) - E{2s, k)cn {2s, k) 


k^ = 


\ - cn {2s, k) 

•Jc]+e,-c, 


(1 - cn {2s, k))^ 

r dt 

E{ip,k) = 

Jo 


1 


9 = 


'2\i 


2 (C2+^2) 

we get the expressions of Bi{s) as following: 


(fl = cos 


-1 


Vl - k^sin^ t 

f5 + 2^C,+f; 


5i( 5) = r I3{t)dt = Qil/i{s) + Dx, 

Jo 

r 7 V 2 

B2 {s) = I 0^{t)dt = -[-35 + (1 - cn (25, k))if/2{s) + i/^ 3 (S)] + D2, 

Jo 


9 


B3{s) 


f 


0^{t)dt 


1 


2g2;t'2 

2 


{k'^cs^ {s, k) + k'\2k^ - l)i/ri(S) - {2k^ -k^- + D 3 , 


B4 {s) = f l2J{t)dt = ^[-^5 + (3 - Ak^)il/2{s) + Ak'\\ + k^Ys - 2k'^E{2~s, k)] + D4, 
Jo 
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where S = ^ ± -^1, F = F((pi,k), D\,D 2 ,D^,D^ are eonstants, and 


Di = -giAi (-), 


£>2 = 
£>3 = 

D,= 


3F , 

(F\ 


— - (1 - c/7 iF,k))ll /2 




9 

2^ - 1)^1 (I) + (2k^ -k^- fc2)e<A.(0 

^f2 r^F /F\ 

^ — - (3 - ( 2 ) - + k^)F + 2k'^E{F, k) 
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